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Abstract 

In this paper we continue previous efforts in the literature to determine phenomeno- 
logical values for the gluon mass by confronting theoretical results obtained in a 
theory of massive gluons with experimental values or results directly referring to the 
nontrivial structure of the Yang-Mills vacuum, e.g. to the presence of the gluon con- 
densate. The decays of heavy quarkonia into 3 gluons and 2 gluons + photon are 
considered in detail as well as the correlators of the electromagnetic current and the 
scalar glueball current. Based on the analysis for the latter quantities a value for the 
gluon mass in the range of 500-600 MeV is estimated from the standard SVZ-value 
of the gluon condensate. 



1 Introduction 



Starting with the original suggestion by Cornwall [[!]], the possibility of a nonzero gluon 
mass has been discussed in various contexts over the past 15 years. It has been argued (c.f. 
0) that the severe infrared singularities of the gluon interaction might resolve themselves 
in the creation of a gluon mass. The emergence of this mass might be closely related to the 
nontrivial nature of the QCD vacuum, in particular to the gluon condensate of Shifman, 
Vainshtein and Zakharov 0. 

However unlike the by now very much established gluon condensate or the constituent 
mass of quarks, the concept of a gluon mass meets a lot of reservation as it touches the 
fundamental principles of gauge invariance and renormalizability. 

As argued by many authors and in particular by ||, the only way to make massive 
nonabelian vector mesons compatible with a decent high energy behaviour, i.e. with renor- 
malizability, seems to be the Higgs mechanism. The effective Lagrangian advocated in 
for a massive gluon theory [] in fact can be understood as a Higgs model with the ram- 
ification that all physical Higgs particles have been removed from the theory by making 
them infinitely heavy. More precisely, this model is the generalization of the SU (2) Higgs- 
model with a frozen radial degree of freedom to the case of the gauge group SU(3). In 
this gauged nonlinear cr-model, the 8 would-be Goldstone bosons provide the longitudinal 
polarization degrees of freedom for the gluons. In the unitary gauge, they are absorbed 
into the redefinition of the gauge field, leaving us with a massive vector boson theory with 
a global SU (3) symmetry. 

Though the conflict between a gluon mass and gauge invariance is avoided in this model, 
the renormalizability has obviously been lost in the way of going from a linear cr-model to 
the nonlinear one. 

One of the main questions then is, how severe the loss of renormalizability really is. 
From the situation of SU(2), where the resulting new divergence is only logarithmic in the 
cut-off or the physical Higgs mass - a fact known as screening theorem of Veltman || -, we 
may hope for a similar behaviour in the SU(3) case. To clarify this point, the systematic 
analysis of |J for the divergencies should be extended from SU{2) to SU(3). 

Another critical issue to the significance of this model as an effective low energy theory 
is the confinement property. In the Higgs phase, i.e. the phase corresponding to super- 
conductivity, magnetic monopoles are confined and the Wilson loop in the fundamental 
representation satisfies a perimeter law. It has been argued in |7|] that at some critical cou- 
pling the model is expected to undergo a transition to a phase with condensed vortices. In 
this phase the 't Hooft loop, i.e. the order parameter related to confinement of monopoles, 
can be shown to satisfy a perimeter law implying - according to the classification of phases 
by 't Hooft - an area law for the Wilson loop in the fundamental representation. 

On the phenomenological side, effects of a finite gluon mass have been considered in 
the past for the standard quarkonium decays into gluons ||, potential models for glue- 

1 note, that we do not follow the suggestion of Cornwall to use the classical equation of motion to express 
the nonlinear er-field in the Lagrangian in terms of the gauge field 
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balls |TIJ, approximations to the QCD-Hamiltonian combined with a variational method 



TT| , a possible relation between the gluon condensate and the gluon mass through QCD 



sum-rules ||12|| , the two-gluon exchange contribution to the Pomeron ||13|| , and the vacuum 
functional of SU(3) Yang-Mills theory 0. 

In a recent analysis ]I5 of J/ip and T decays it was investigated in detail whether 
the deficiencies between experiment and theory for the total decay rate and the inclusive 
photon spectrum can be understood as being due to relativistic corrections or alternatively 
due to a finite gluon mass. The conclusion of [ 15 1 is that the relativistic corrections needed 
to bring experiment and theory into agreement do not scale in the expected manner when 
going from c-quarks to b-quarks. The alternative with a gluon mass leads to values of 
.QQGeV and l.llGeV for J/ip and T decays, respectively. The authors interpret the rise 
of the gluon mass parameter as an indication that without phase-space limitation the true 
gluon mass would be of the order of 1.5GeV and is reduced to effectively smaller values 
due to phase space limitations in T and J/ip decays. 

As a continuation of previous efforts to determine phenomenological values for the gluon 
mass, we treat in this paper the following three subjects. First we complete in section 2 the 
analysis of |Tj| by taking into account a finite gluon mass not only in the phase space but 
also in the amplitudes for J / ip and T decays. For this purpose we derive the analytical form 
of the Dalitz distributions for {QQ)\ — > 3g, (QQ)\ — > 7 + 2g and (QQ)l — > 2g for massive 
gluons to leading order in a s , thus extending the classical formulas originally derived for 
positronium decays (e.g. ref. (16|]). Using these results we determine by a comparison with 
the results of |L5| corrected values for the gluon mass. 

As a further basic process affected by a finite gluon mass we consider in section 3 the 
e + e~ annihilation into mass-less quarks. The analytic form of the vacuum polarization 
to order a s with a massive vector boson has been known for a long time |]17|| . Analyzing 
the dispersion relation for the derivative of the vacuum polarisation function, a systematic 
expansion in powers of M 2 /Q 2 can be derived for the Euclidean region. This expansion 
can serve as a basis to relate the gluon mass M to the gluon condensate by equating the 
1/Q A term in this expansion with the 1/Q A term in the operator product expansion as used 
before in ref. jjEJ. As it turns out, the expression obtained in this way is different from 
the expression in ref. [|12| , where a Mellin-transform technique has been used to isolate the 
1/Q A term in the loop diagrams. Our relation can be used to derive a simple expression 
for the gluon mass in terms of the gluon condensate. 

In a similar spirit we deal in section 4 with the two-point function for the scalar glue-ball 
current instead of the electromagnetic current. Matching the result obtained for this two- 
point function in the massive theory with the result obtained from the operator expansion 
in the massless theory, a similar, albeit not exactly identical, relation between the gluon 
mass and the gluon condensate is obtained. 

In section 5 finally the results are summarized and a few remarks concerning possible 
future directions of this research are made. 
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2 Heavy quarkonia decays into gluons 



In this section we present our results for the decay rates of triplet (QQ)\ and singlet (QQ)l 
S-w&ve quarkonia states into massive gluons. The calculation uses the standard approx- 
imation of neglecting the dependence on the momenta of the heavy quarks. Accordingly 
we can write: 

dT((QQ) - ..) ~ |^(0)| 2 rfr(Q + Q - ..) (1) 

where "0(0) represents the spatial QQ wave-function at the origin and dT(Q + Q — > ..) is 
calculated with Q and Q at rest. 

We use 7] = (M/2m) 2 to parametrize the dependence on the ratio between gluon mass 
M and quark mass m. On the level of our approximation, 2m of course can be identified 
with the mass of the respective singlet or triplet quarkonium state. 

We start with the Sq — * 2g decay rate which is given by: 

r = r ((l- 4^/(1 -2r ? ) 2 )(l-4r / ) 1 / 2 (2) 

where T is the decay rate for M = and the gluon mass effects on the amplitude 
squared and the phase space have been separated for convenience. As can be seen from 
eq.(||) the amplitude squared decreases with increasing M, though compared to the phase 
space the effect only starts at one order higher in rj . 

For the Dalitz distribution of 3 Si — > ggg we obtain the basic formula: 



dT( 3 S 1 ^ggg) = (3) 
r 5 (—2 (1 — 3 rj) + x\ + X2 + £3) dx\ dx^ dx-$ , 

( — 9 + 7T 2 ) Xi 2 X2 2 X3 2 



3 v " \ 8 4 J 
+ (xi 2 + x 2 2 + x 3 2 ) (l - 14 rj + 48 rf + 25 rj'" 

-2 (x! 3 + x 2 3 + x 3 3 ) (l - ^ - 3 rj 2 



+ (x x A + x 2 4 + x 3 4 ) (l + l 
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where the variables X{ are related to the gluon energies E± by Xi = Ei/m — 2rj. Similar to 
before T represents the integrated decay rate for mass-less gluons. 

Apart from the constant term proportional to rj , eq. (Q) has a similar form as in the mass- 
less case, i.e. the numerator of the distribution can be written as (x 2 P(xi) + cycl), where 
P(x) is a second order polynomial in x. For rj = 0, P reduces to the form P(x) = (1 — x) 2 , 
well known from the classic positronium formula. 
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When one of the gluons, say gluon 1, is replaced by a photon, the corresponding ex- 
pression reads with z = E y /m replacing x±: 



dT( 3 S 1 ^ 1 gg) = 

T 5 (—2 (1 - 2 77) + z + x 2 + x 3 ) dzi dx 2 dx 3 



(4) 



(-9 + 7T 2 ) 

4"( l - 2J r) f 1 - 2 ") 2 

+^32 (1-2,7) (l-£) rf 
+z 2 (I-277) (l -6 77 — 6 rf) 



z 2 x 2 2 x 3 2 



-2z A 



+ (x 2 2 + a;3 2 ) (l -8t7 + 22t7 2 + 8?7 
-2 (x 2 3 + x 3 3 ) (l 



+ (s 2 4 + s 3 4 ) + ? 



(x 2 2 x 3 2 ^ 77] 
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Figure 1: 3 Si — > (7(7(7 decay rate as a function of rj = M 2 /4m 2 , normalized to the decay 
rate for mass-less gluons. The dashed line represents the case when the gluon mass is only 
taken into account in the phase space. The dotted lines indicate how the gluon mass has 
to be changed to reach the same suppression when taking the gluon mass into account also 
in the amplitude 




Figure 2: The analogue of Fig. 1 for 3 Si — > 7 + gg 
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Using these expressions and integrating over phase space we can now make a simple 
analysis to estimate new values for the gluon mass which take into account both phase- 
space and amplitude effects. The results are shown in Figs. 1,2, where the integrated decay 
rate is plotted as a function of M 2 /4m 2 . As can be seen, the inclusion of the gluon mass 
in the amplitude leads to an increase of the decay rate. This feature might be expected on 
the basis that for massive gluons there are in addition the channels involving longitudinal 
gluons. Using the values for the gluon mass from ref.(||15||) i.e. M = .66GeV and M = 
1.17GeV for J/ip and T, respectively, we conclude that in order to reach the same level of 
suppression the gluon mass values change to M = .76(.79)GeV and M = lAl(lA)GeV, 
where the numbers in brackets refer to the case where the 7 + 2g reaction is used to make 
the matching. 

Another quantity of interest is the photon energy spectrum. The analytic form is 
given in the Appendix. Examples for the effect of the gluon mass on the shape of this 
distribution are shown in Fig. 3. As can be seen, the larger gluon mass needed to produce 
the same amount of suppression when the gluon mass effect is taken into account in both 
the amplitude and the phase space, leads to a sharper cut-off at the end of the spectrum. 
However, in this region we cannot take the massive gluon model really serious, and there 
are other effects like the Fermi-motion of the quarks to which this behaviour is sensitive 
to. 



3 e + e annihilation 

As a further basic process we consider now the leading order effect of a massive gluon on 
the electro-magnetic vacuum polarisation, i.e. the diagrams given in Fig. 4. To simplify 
matters, we take the massless quark limit. A systematic expansion in terms of M 2 /Q 2 
can most easily be derived starting from the dispersion relation for the derivative of the 
vacuum polarisation function. 



D(Q 2 ) = Q 2 ds-^- 2 \ (5) 




The imaginary part R for the diagrams in Fig. 4 has been known for a long time |17| , and 
for x = M 2 /s < 1 is given by: 
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Figure 3: The normalized photon energy spectrum for T — > 7 + gg as a function of 
z = Ej/m. The pair of gluon mass values chosen, i.e. M = lAlGeV and M = l.l7GeV, 
are the ones determined above to reproduce the experimental suppression of the total 
decay rate. As before, the dashed line is the curve involving the phase-space effect only. 
For comparison also the spectrum for zero gluon mass is shown (dotted line) 




Figure 4: The leading order diagrams with a gluon of mass M 
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R{x) = (6) 
N 2 — 1 3 

° , r ( 2x-5x 2 + 2x ln(x) + 3x 2 ln(x) 

AN C 2 

-2 (1 + xf (ln(x) ln(l + x) + Li 2 (-x)) ) 

For convenience the normalisation has been factored out, i.e. N c a s /ir for a single quark 
flavor with unit charge and N c colours. Eq.(^) also determines R for x > 1 through the 
symmetry relation: 

R(x) = ( 1 - x 2 ) R(0) + x 2 R(-) (7) 
The function i?(x) has the small- x expansion: 



R(x) I (^2_i) \ = R(x) = l+(Y / a k x k ) + (Y / hx k ln( 



fc=2 / \fe=3 



As a consequence, the following expansion for D(Q 2 ) in terms of M 2 /Q 2 can be derived: 

DW .),(^)i= ( 9 ) 

1 + ^2 C/'d..^-^ 



<5 2 \-A) x 2 / 

( M 2 V ( ,„,,M 2 , n ri £(x)-l n /"i i?(x)-l-a 2 x 2 
+ 7^ a 2 (2 ln(— ) + 1) - 2 / dx 2 / dx ■ 



Q 2 J \ Q 2 Jo x Jo x 3 J 

M 2 3 

+0(^) (10) 



Q 

Numerically we have: 



a 2 = -1 

/ 1 (fa M^I = _ 7+ il! + 2C (3) = -0.209 
Jo x 9 

i (fa fl(x)-l = -16 +4 c(3) = _ 0-525 

n X 2 3 

^(x)-l + x 2 57- 8vr 2 + 36C(3) 

dx = — = 1.184 

x A 18 



Equation (|S]) should be compared with the result from the operator-product expansion: 
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The first term represents the perturbative contribution from the diagrams in Fig.4 with 
massless gluons, and the second term is the contribution from the gluon condensate as 
given in ref. || (in our normalization). Equating the 1/Q A terms in (|9|) and (0), we 
obtain: 



< — GG >= — M 4 ( ^ - 8 C (3) - 2 In, ^ ) ) , 1 3 ) 

7T 7T IDF 




The corresponding expression derived in |T^] reads: 



— GG >= — 3 ^ ~ a M 4 (— + — C (3) -*J£-JL H *)-* H "L* 
n Ti 16 7i 2 I 8 3 sv ' 27 18 K Q 2 ' 9 K Q 2 ' 



which has no resemblance to what we have obtained. 

At present the origin of this discrepancy is not satisfactorily explained. However, since 
our own calculation only uses standard Feynman diagram techniques and has been checked 
by an independent calculation where the diagrams in Fig.4 are calculated by a numerical 
procedure along the lines of ref. [TS| , we suspect that the Mellin transform technique used 



in ref. 0] has not been properly applied. 

A Q 2 -independent relation between M and the gluon condensate can be obtained from 
eq.(|T3"D by using for a s /7T the standard expression 4/ /3oln(Q 2 / A 2 ) for the running coupling 



with p = N c ll/3 and evaluating eq.(|13|) to leading order in ln(Q ). The result is: 



M 4 = — 7r 2 - — ^ — - < — GG > (15) 
9 (A 2 -l) vr 1 J 

yielding M ~ 570MeV for N c = 3 and the standard SVZ-value for the gluon condensate 
< s& GG >= (.33Ge^) 4 . 



4 Scalar glue-ball current 

Instead of the electromagnetic current we consider in this section the scalar glue-ball current 
defined in terms of the gluon field-strength tensor G a by 

] = G^ u G afW (16) 

For ease of writing we use j to denote the current with a unit normalization factor and 
j for the current with a normalization factor a s /iT which turns j into a renormalization 
group invariant quantity f\. 

2 For a leading logarithm analysis it will be sufficient to use a s /ir instead of the more generally valid 
factor 40(g)/ (fog) 
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As a first step we calculate the diagram in Fig. 5 for j in the massive theory. Defining 
fr GG (Q 2 ) for euclidean q 2 = —Q 2 by: 



N 2 - 1 

<jj>=i-^*GG(Q 2 ) (17) 



we obtain: 



^,m=(\{Q 2 + ^M 2 ) 2 + M^ ( f'da Hl + a( \ ^ Q2 )) (« 



2^ l ) \Jo M 2 

where = means equality up to polynomial terms which do not survive a three-fold subtrac- 
tion needed to remove the UV-divergencies. From (|T^) we obtain the following expansion 
in terms of M 2 /Q 2 : 

,1 m 2 M 2 M 2 

QH ^ naa ^. 3 ^ + 12( ^ r + 0{ p (19) 

which is free of logarithmic factors ln(M 2 /Q 2 ) up to order l/Q 4 . As visible from this 
expression, the normalization of 7t GG has been chosen such that the third derivative equals 
to one in the massless limit. 



M 




2/ii - k l k 2 9^i^ 2 ) 



Figure 5: The leading order diagram for the scalar glueball current with a gluon of mass 
M, and the expression for the vertex. 

As a second step we make a renormalizat ion-group analysis for the coefficient function 
of the gluon condensate < j > in the operator- product expansion of the two-point function 
n GG for the scalar current j = (a s /n) j: 

kgg = C pert (Q 2 ) 1 + C GG (Q 2 ) < j > + ... (20) 
The coefficient function C GG has the following expansion in a s /n : 
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Caa{Q 2 ) = ^ Co + (^) 2 Cx(^) + ... (21) 



7T 7T Q z 

where Co is a constant whose value can be found from applying Wick's theorem to the 
product of currents j(x)j(0): 



Independence of Cqg from the renormalization point \i requires the coefficient C\ in ( pT|) 
to be a polynomial of degree i in the variable In^/Q 2 ). Choosing /i 2 = Q 2 and using the 
1-loop expression for a s (Q 2 )/ir, we obtain: 

where = H/3iV c 

The desired relation between M and the gluon condensate is now obtained by matching 
the 1/Q 4 terms in eq.(|23|) and (a s (Q 2 )/n) 2 x eq.tfHty. The result is: 



Comparing with the result obtained in section 3 in eq. (|i5|) , we recognize that we have 
almost obtained the same result except that the numerical coefficient is smaller by a factor 
of 2. Thus the previous estimate for M is reduced from 570MeV to M ~ A80MeV 

It is interesting to compare eq.(24) with the result of ref.( ||14|| ) obtained by a variational 



investigation of an ansatz for the vacuum functional. Again both results have the same 
structure and color dependence (large N c ). However, here there is a large difference in the 
numerical constant, i.e. 120 for ref.(|l4|) instead of 11/9. As a consequence their mass 
estimate is larger by a factor ~ 3. 
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5 Summary and Outlook 



In this paper we have continued previous attempts to estimate values for the gluon mass 
by confronting results obtained in leading order calculations for massive gluons with ex- 
perimental results; either by a direct comparison as for J/ip and T decays or in a more 
indirect way by matching against nonperturbative contributions in the correlator of the 
electromagnetic current or the scalar glue-ball current. For the decays of heavy quarkonia 
the Dalitz distributions involving 2— or 3— massive gluons were derived, extending thus 
the basic formulae originally derived for positronium to this more general case. We then 
analyzed the leading QCD correction of the electromagnetic vacuum polarisation and de- 
termined the coefficients in an expansion in 1/Q 2 which are relevant for a comparison with 
the leading nonperturbative contribution for this quantity. We found a different result 
for the l/Q 4 -coefficient as previously given in the literature and claim that our result is 
correct as it has been checked by a completely different method than the one presented in 
this paper. Our result was then used to derive a simple renormalization-group invariant 
expression between the gluon mass and the gluon condensate. Finally a similar kind of 
analysis was applied to the correlator of the scalar glue-ball current where gluon mass 
effects already appear on the 1-loop level. Again this leads to a relation between the gluon 
mass and the gluon condensate which is renormalization group invariant although with a 
slightly different numerical coefficient. 

Taken at face value, there is a large spread between the mass estimates from heavy 
quarkonia decays, in particular the T-decays where we can expect gluon mass effects to 
have fully developed, and the values obtained from the analysis of the correlators. Certainly 
the value of M ~ lAGeV obtained from T-decays should be considered as an upper bound, 
since the starting values of |15j for the gluon mass have been obtained under the assumption 
that the main suppression comes from the gluon mass, not the relativistic corrections. As 
argued for instance in [[L9| , the pattern of suppressions in the total J / ip and T decay rates 
can quantitatively be reproduced if the momentum dependence of the quarks is taken into 
account. Whether this straightforward incorporation of relativistic effects is legitimate is, 



however, still under debate (see [fig]), and therefore a substantial gluon mass - though 
smaller in magnitude than the values discussed here - may indeed be required by the data. 
The results from matching the leading nonperturbative correction in the correlator for the 
electromagnetic and the scalar glueball current are equal to each other up to a factor of 2 
in the relation for M 4 . Although this has little effect on the numerical value of the gluon 
mass M, the dependence on the process remains to be explained. In this context it would 
also be interesting to derive the expansion in powers of M/Mqq of the total decay rates 
of heavy quarkonia discussed in section 2 and match them likewise against the leading 
nonperturbative correction from the gluon condensate. 
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A more complete analysis should also cover the 1/Q 2 corrections obtained in the massive 
theory, which seem to have no correspondence in the operator-product expansion. It has 
been argued in ref.(||12||) that these corrections cancel against corrections from the formation 
of a vortex condensate. We believe that this proposition requires further examination. 
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7 Appendix 

For the distribution of the photon in 3 Si — > 7 + gg the following expression is obtained 
from eq.(||) by integrating over the gluons. 



dT 

dz 



oc 



X-L- — X- 



(25) 



+ 



ln{x + / xJ) 



[8 (—1 + 2/7) 2 (2-4// + 7; 



z 2 (-2 + 4 77 + zf 
+8 (-1 + 2 77) (5 - 12 77 + 10 rj 2 + 2 r/ 3 ) z 

+2 (-1 + 2 77) (-17 + 10 7] + 6 rj 2 ) z 2 

+2 (-5 + 277 + 2r7 2 ) z 3 } 

z 2 {—2 + 4ri + z) v 
+4 (-1 + 27?) (3-4r/ + 2r? 2 + 2r/ 3 ) z 
+2 (7 — 18 77 + 10 r? 2 + 10 ry 3 ) z 2 
+4 (2 + 7/) (-I + 277) ^ 3 + (2 + 77) z 4 ] 



Here z = E^/m in terms of the photon energy E^, and 
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delineate the boundary of the Dalitz plot in the plane (X2, z) with X2 = E^/m — 2r). 
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